Spatially Resolved Star Formation History Along the Disk of 
M82 Using Multi-Band Photometric Data 

L. H. Rodriguez- Merino, D. Rosa-Gonzalez & Y. D. Mayya 
INAOE, Luis Enrique Erro 1, Tonantzintla, Puebla, Mexico, CP. 72840. 

linoOinaoep .mx 



ABSTRACT 

We present the results on the star formation history and extinction in the disk of M82 over 
spatial scales of 10" (~180 pc). Multi-band photometric data covering from the far ultraviolet 
to the near infrared bands were fitted to a grid of synthetic spectral energy distributions. We 
obtained distribution functions of age and extinction for each of the 117 apertures analyzed, 
taking into account observational errors through Monte-Carlo simulations. These distribution 
functions were fitted with gaussian functions to obtain the mean ages and extinctions along with 
errors on them. The analyzed zones include the high surface brightness complexes defined by 
lO'Connell fc Mangano ( 19781) . We found that these complexes share the same star formation 
history and extinction as the field stellar populations in the disk. There is an indication that the 
stellar populations are marginally older at the outer disk (450 Myr at ~3 kpc) as compared to the 
inner disk (100 Myr at 0.5 kpc). For the nuclear regions (radius less than 500 pc), we obtained 
an age of less than 10 Myr. The results obtained in this work are consistent with the idea that 
the 0.5-3 kpc part of the disk of M82 formed around 90% of the stellar mass in a star-forming 
episode that started around 450 Myr ago lasting for about 350 Myr. We found that field stars 
are the major contributors to the flux over the spatial scales analyzed in this study, with stellar 
cluster contribution being 7% in the nucleus and 0.7% in the disk. 

Subject headings: galaxies: individual (M82)- galaxies: formation- galaxies: star clusters 



1. Introduction 

The M82 galaxy (NGC 3034) is one of the 
most studied extra-galactic objects in the local 
universe. It is an edge-on spiral galaxy clas- 
sified recently with a morphological type SBc 



( Mavva et al. 20051) . The sma ll distance from 



the M ilky Way to M82, 3.63 Mpc (|Freedman et al 



1994I ). has allowed to obtain images with high spa- 
tial resolution, which were e mployed to study its 
strong star-forming activity (Solinger et al. Il977t 



% ■ < , , 

Telesco et al. Hl99lHForster Schreiber et al. 120011) . 

its compact star clusters ( CSCs) distributed along 
the disk of the galaxy (lO'Connell fc Mangano 



19781: lO'Connell et al. lll995Hde Griis et all 12001 



Melo et al.l2005HSmith. et al.l2006ll2007t lMavva et al. 



20081 : iKonstantopoulos et al.l 120091) . its biconi- 
cal outflow which produces a filamentary struc- 



ture of sevexal^fl^ars^c4ong_(Lvnds & Sandagc 



19631: ICappi et al.l 1 19991: iHoopes et al. I 120051) . 



and the properties of the neutral and molecu- 
lar gas (IShen fc Lo Ill995t iNeininger et al. |[i~99j 



Matsushita et al. I 120051 ) . Particularly, the ob- 



served bridge of neutral hydrogen, which is con- 
necting M82 with M81, has been used as a 
test of the interaction between t hese objects 



(ICottrell I Il977t lYun et al. I Il993l: ISofud Il998l: 



Chvnoweth et al."T2008 ). There were attempts 



to determine the exact epoch of interaction, with 
the estimated epoch dependin g highly on t he as- 



Yun et al. I (Il993l) 



sumptions in the simulations, 
reproduced the observed HI plumes around M82 
by a model where the last passage o f the compan - 



Sofuel (|l998h 



ion occurred around 300 Myr ago. 
estimated the epoch of interaction as around 1 Gyr 
ago. This was obtained assuming that the funda- 
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mental stellar disk of M82 was blown away during 
the interaction, and that is the reason for the ob- 
served nearly Keplerian rotation curve. Given 
that t he disk is rich of a population of AGB 
stars ( Davidge 20081) — the products of post- 
interaction star-formation episode — the timescale 
seems to be longe r than the 300 Myr estimated by 
Yun et al. I (|l993h . We hence use an intermediate 
age of 500 Myr as the epoch of interaction, for the 
sake of discussion in this work. 

O'Connell fe Mangano (Il978h using the best 



of fos sil starburst. More recently, ISmith et al.l 



available photographic images of M82 at that time 
analyzed the brightest areas of M82. They identi- 
fied 8 complexes, which they denoted by the let- 
ters A to H, and pointed out that most of these 
complexes contain knotty structures. They found 
that along the disk (r >1 kpc) there is no trace 
of complexes with on-going star formation. They 
also noted that the nuclear part of the galaxy (ra- 
dius less than 500 pc) presents a strong, recent 
star formation acti yity, with th e complex A being 

" (11980 . Il993h 



the most intense. Riekeetal 



ana- 



lyzed the multi-band properties of this zone, and 
determined a n age of less than 30 Myr for th e nu- 



clear region. iForster Schreiber et al. ( 200l[ ) car 



ried out a spatially resolved study of stellar ages in 
this region, concluding that the nuclear zone had 
more than one burst of star formation in the last 
10 Myr. These authors have noted that the patchy 
appearance of the nuclear zone is due to dust ob- 
scuration with more than 4 mag of visual extinc- 
tion if the majority of the obscuring dust resides in 
a foreground screen and 43-52 mag if dust is mixed 
with the stars. Sub-arcsec spatial resolution im- 
ages off ered by the Hubbl e Spac e Telescope (HST) 
allowed O'Connell et al. (1995) to discover more 
than a hun dred compact stellar clusters in the nu- 
clear zone. ISmith et al.1 (|2006h . using HST/STIS 
spectroscopy, found that these clusters are young 
and are formed as part of the nuclear starburst 
that have been taking place over the last 10 Myr. 

The complexes outside the nuclear zone are 
not as extensi vely studied as those in the star- 
burst zone. Ide Grijs et al ] (|200lh . using the 
HST/WFPC2 images, found plenty of compact 
stellar clusters in the complex B, all older than 
those in the nuclear region, but younger than 
around 500 Myr. They found that the star forma- 
tion rate in complex B was as high as the nuclear 
complex at the present day, thus giving it the name 



(|2007l ) has derived a mean age of around 150 Myr 
for 35 clusters in the complex B. Complex F is 
the only other zone wh ere studies in some detail 



have been carried out (Smith fc Gallagher 2001 



Bastian et al.l 120071) . These studies gave an age 



of less than 100 Myr for this complex, with no 
evidence for recent star forma tion. More recently, 
iKonstantopoulos et al. (|2009h have derived spec- 
troscopic ages of around 50 compact star clusters 
distributed over the disk, with the median age of 
the clusters being 150 Myr. 

In order to explain the presence of patchy 
bright complexes in the nucl eus and the disk, 



O'Connell fc Mangano (|l978l ) invoked a model 
where the star formation was restricted to only 
the bright complexes. They suggested that after 
the tidal encounter of M82 with the galaxy M81, 
some interstellar material was stripped from M82. 
When this material fell back into the galaxy, it 
produced the classical bright complexes with en- 
hanced episodes of star formation. If this model 
is correct then the main population of stars in the 
disk surrounding the bright zones should be older 
than the bright zones themselves. 



On the other hand, iMavva et al. J (120061) using 
photometric, dynamical and chemical properties 
of the disk between 1 and 2.7 kpc radius proposed 
a violent star-formation episode throughout the 
disk of M82 that started around 0.8 Gyr ago, and 
lasted for a few hundred million years. This 
model of star formation reproduces well the rel- 
ative strengths of age-sensitive spectroscopic fea- 
tures in the blue part of long-slit spectra along 
the major axis, as well as the observed value of 
M/Lfc, metallicity and enrichment of a-elements. 
They also noted that the galaxy does not possess 
a massive old underlying stellar disk, characteris- 
tic of disk galaxies. More than 90% of the stel- 
lar mass in the disk was formed in this episode, 
implying that the star-formation episode is not 
jus t restricted to a few patches a s was suggested 
by lO'Connell fc Mangano ( 1978 ). instead it was 
disk-wide. The b asic tenets of thi s scenario were 
recently tested by iDavidge" (|2008f ) using near in- 
frared color-magnitude diagrams (CMDs) of the 
resolved populations in the disk. He detected 
asymptotic giant branch stars (AGBs) of age of 
a few hundred million years which are distributed 
uniformly in the disk, up to a radial distance as 
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large as 12 kpc, suggesting that the star forma- 
tion episode pervaded the entire disk. More im- 
portantly, red super giants (RSGs) were absent in 
the disk outside a radius of 0.5 kpc, confirming 
the cessation of star-formation episode ~ 100 Myr 
ago. 

While the bursty nature of star-formation his- 
tory in the disk of M82 is by now well established, 
there is still room for improvement in the tim- 
ing of the commencement of the disk-wide star - 
formation. In the study of iMavva et al. I (|2006l ). 



the uncertainty in age estimation comes from the 
errors involved in determining the stellar mass of 
the disk (to calculate the quantity M/Lfc), which 
was obtained using a dynamical model to fit the 
observed rotation curve of the disk. Given that 
M82 has suffered a major change dynamically as a 
result of its interaction with the members of M81 
group, and the evidence for the existence of gas 
above the plane, some of the intrinsic assumptions 
made to determine the mass, such as planar cir- 
cular motions, may not be completely valid. An 
underestimation of the mass would push the com- 
mencement age to be as old as 1 Gyr, whereas 
if the dynamical masses were over-estimated, the 
mean ages could be younger than 0.5 Gyr. On the 
one h and, the analysi s of the AGB stars carried 
out bv lDavidge" ( 2008h promises to provide an ac- 
curate clock for age dating. However at present its 
utility is limited because of the poor knowledge of 
this evolutionary p hase. Most recent models of 
iMarigo et all (|2008l ) indicate age >200-300 Myr, 
with peak of the AGB population occurring at age 
>500 Myr. Hence alternative techniques are nec- 
essary to determine the age of the stellar disk in 
M82. 

In recent years, the SEDs formed from pho- 
tometric data covering a wide wavelength in- 
terval are increasingly used in determining the 
age, extinctio n and metallicity of stellar systems. 
iBianchi et al. (2005) using ultraviolet (UV) im- 
ages from Galaxy Evolution Explorer (GALEX) 
and optical images from the Sloan Digital Sky Sur- 
vey (SDSS), derived ages and reddening (among 
other parameters) of 160 pc scale stellar complexes 
in M101 and M51. iKaviraj et "all (|2007f ) estimated 
age and metallicity of 42 low-reddening globular 
clusters in M31 using the p hotometry in the FUV , 



NUV, U, B, V, R, I bands. iBridzius et aTl (|2008h 



explored the capability of the UBVRIJHK photo- 



metric system to estimate the age, metallicity and 
color excess of star clusters. These studies have 
illustrated that the long base line in wavelength 
helps to disentangle the effects of reddening from 
the evolution. SEDs can also be used to get a 
measure of the metallicities in old, low-reddening 
systems such as globular clusters. However, the 
photometrically-derived metallicities do not reach 
the accuracies that are achievable using spectro- 
scopic data. 

In the present study, we exploit the potential of 
present-day SED-fitting techniques to determine 
the ages of the dominant stellar populations in 
the disk of M82, and in the process spatially map 
the star formation history. The available photo- 
metric data allow us to analyze the star forma- 
tion history over spatial scales of ~ 180 pc, the 
size limitation basically coming from the 5" beam 
of the GALEX images and the requirement that 
measured fluxes in each aperture have the enough 
signal to noise (S/N) ratio. In e arlier studies 
(jMavva et al. 1 120061; [Davidge II2008I ). star forma- 
tion history was inferred by analyzing data over 
spatial scales of 1-2 kpc, so the present study rep- 
resents an improvement by a factor of around 10 
in spati al resolution. Typically, the com plexes de- 
fined bv lO'Connell fc Manganol (1 19781) have spa- 
tial extensions comparable to our resolution ele- 
ment, with the exception of complex B, which is 
several times larger than our aperture. The use 
of imaging data to determine ages also allows a 
comparison of star formation history of the bright 
complexes with those of the surrounding disk. 

In §2, we describe the observational data used 
in this work. In §3, we describe the method used 
to infer the age of the stellar population and ex- 
tinction. In §4, we describe the results obtained 
for the named complexes. In §5, we analyze the 
spatial distribution of the stellar population age 
to infer the star formation history in the disk and 
nucleus of M82. In §6, we determine the cluster 
formation efficiency and in §7, we discuss our re- 
sults. Concluding remarks are given in §8. 

2. Observed Spectral Energy Distribu- 
tions 

The most important part for deriving the spa- 
tially resolved star formation history from photo- 
metric data is the availability of high S/N ratio 
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Fig. 1. — F555W image of M82 from HST. The positions of the 117 apertures (5" radius) defined in this 
work are marked with circles. The classical M82 complexes A, C, D, E, F, G, and H correspond to the 
apertures 57, 65, 59, 64, 72, 80 and 52 respectively. The complex B contains the apertures 27, 29, 32, 35, 
40, 41, 45 and 50. 



images covering the wavelength interval from the 
far ultraviolet (FUV) to the near infrared (NIR). 
The S/N ratio in the analyzed zone should be bet- 
ter than 20 in order to obtain reliable ages. 

The UV part of the SEDs is based on the 
GAL EX far-UV (FUV) an d near-UV (NUV) im- 
ages ( Morrissev et al.ll2007 ) . The optical images in 
the u, g, r, i and z bands were retrieved from the 
SDSS database (jAbazaiian et al.l 120031) . For the 



near infrared data, we used the J, H and K images 
obtained with the Cananea Near-Infrared Camera 
(CANICA) attached to the 2.1 m telescope at Ob- 
servatorio Astrondm ico Guillermo Haro (OAGH, 
iMavva et ai~ll2005l) . 

Among all these dataset, the GALEX images 
have the poorest spatial resolution with a point 
spread function of full width half maximum of 
5". These images set the minimum spatial scale 
of our analysis, which corresponds to ^100 pc in 



M82. We carried out photometry in fixed aper- 
tures of 5" radius. The apertures were placed 
so as to sample the bright co mplexes defined by 
O'Connell fc Mangano I dl978h . as well as the dif- 



fuse part of the disk. Figure Q] displays the F555W 
image of M82 from HST, showing the positions 
of the 117 selected apertures. In Table [TJ we 
show the association of o ur ap erture numbers and 



O'Connell fc Mangano I (jl978l) complexes. Aper- 
ture numbers increase from east to west along the 
major axis of the galaxy. Apertures 50 to 80 (see 
Figure Q} lie within the 500 pc of the center of the 
M82 and belong to the nuclear starburst. 

We aimed that the ages inferred from the SED- 
fitting technique are good to 15%. We found that 
the photometric accuracy of each band used for 
the construction of the SED should be better than 
5% in order to achieve this. This corresponds 
to limiting surface brightnesses (for a S/N ratio 
of 20), averaged over an aperture of 5" radius, 
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Fig. 2.— Mean SEDs across the M82 disk. Solid 
line represents the mean of SEDs for selected nu- 
clear apertures. Dashed and dot-dashed lines are 
the mean SEDs for apertures located in the north- 
east and south-west of M82 respectively. The 
mean SEDs were normalized in the K-band flux 
and the vertical bars represent the flux dispersion 
(la). 



of: FUV=28.01, NUV 
r=18.66, i= 18.23, z= 



26.85, u=19.34, g=19.97, 
16.69 (all in units of AB 
magarcsec- 2 ), J=23.16 (Vega), H=19.89 (Vega) 
and K=19.23 (Vega). M82 is brighter than these 
limits up to a galactocentric radius of 3 kpc in all 
bands, except for 8 apertures located in dusty fea- 
tures, whose FUV and u magnitudes were fainter 
than the limiting magnitudes. We eliminated 
these apertures from rest of the analysis. Once the 
apertures were defined, photometry is carried out 
by summing all the sky-subtracted flux inside the 
apertures. The sky value for an image is obtained 
as the mean value of several selected areas outside 
the main disk of M82. Obvious sources, such as 
stars and background objects were avoided in the 
definition of sky zones. Note that we subtract sky, 
not the background, thus ensuring that our pho- 
tometry measures the flux emitted by all the stars 
in the disk. 

The error (A/) in the flux (I) within the aper- 
ture was calculated by using, 



AT 
T 



counts+sky 
gain 



Npix * <7 2 



BG 



(i) 



counts 

where, Npix is the number of pixels inside the 



aperture, counts and sky are the total number of 
counts in selected areas in the disk, and in ar- 
eas free of sources, respectively. The obg in the 
equation is the RMS error per pixel, which was es- 
timated using the sky regions in each image, and 
gain is the number of electrons per data counts 
which is provided by the GALEX and C A NIC A 
web pages, and in the case of SDSS, it is found in 
the related fpAtlas file. The instrumental errors in 
the estimated aperture fluxes calculated using the 
above formula is less than 1% in all bands except 
in the bluest bands (FUV, NUV and u). In the 
FUV -hand the errors for majority of the apertures 
lie between 5-10%,while errors lie between 1-5% 
for majority of the apertures in the NUV and u 
bands. On the other hand, external errors, which 
are mainly determined by the errors in flux cal- 
ibration between the three datasets, could be as 
large as 5%. Both errors were added in quadra- 
ture to estimate the total error on each datapoint. 

In each image, we obtained the sky-subtracted 
fluxes in the selected apertures to construct 117 
spectral energy distributions (SEDs). Figure [5] 
shows the mean of nuclear SEDs (solid line), 
north-east SEDs (dashed line), and south-west 
SEDs (dot-dashed line). The SEDs were normal- 
ized to the K-band flux, to facilitate the compar- 
ison. The dispersion of individual SEDs from the 
mean is shown by the vertical bars at each filter 
band. It can be seen that the SEDs of north-east 
and south-west follow a similar trend with the lat- 
ter being marginally bluer. The behavior of the 
mean SED of the central apertures is different. It 
is redder and peaks at longer wavelengths as com- 
pared to the apertures in the disk. Understanding 
these differences in terms of different star form- 
ing history and extinctions, is one of the topics of 
discussion in the current paper. 

3. The SED Fitting Method 

The method we followed to determine the age 
and extinction of the stellar population has six 
steps: 

1. Construction of synthetic SEDs: 

The first step in obtaining the age and 
extinction of the stellar population is the 
construction of synthetic SEDs, suitable 
for comparison with the observed SEDs. 
We used a set of 99 simple stellar popula- 
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Fig. 3. — Age and color excess obtained for complex B (aperture 27). The Monte-Carlo simulation provides 
several best-fit models. The fitted gaussian (dashed line) gives us the value of the age (left) and extinction 
(right) and their corresponding errors. 



tion (SS P) models computed by the Padova 
Group (|Bressan et all Il993l ) 0. The set of 
models cover evenly - in logarithmic scale - 
an age range from 1 Myr to 10 Gyr. These 
models use S alpeter's IMF a n d hav e solar 



metallicity. iMcLeod et al. I (|19931 ) using 



nebular lines found that the metallicities of 
the nuclear regions of M82 are close to solar. 
There arc no metallicity measurements in 
the disk regions due to the lack of emission 
lines there. Following the observed trend 
of metallicities staying constant or decreas- 
ing with galactocentric distance in late-type 
galaxies, we expect solar or sub-solar metal- 
licities in the disk. However, we used solar 
metallicity models for all the regions. 

In order to compare synthetic SSP spectra 
with the observed SEDs, we construct a grid 
of synthetic colors by integrating each syn- 
thetic spectrum over the response curve of 
the filters used in the observations described 
in the previous section. 



most recently updated models kindly provided by A. Bres- 
san. 



Reddening the synthetic SEDs: 

The synthetic S EDs w ere reddened using the 

Calzetti et al.l ( 19941 ) extinction law. We 



also explored the result s obtained if we mak e 
use of the Cardelli law ( Cardelli et al.|[l989l) . 
Once the relative reddening is fixed by as- 
suming an extinction curve, the total extinc- 
tion is computed by varying the color excess, 
E(B-V), from 0.0 to 2.0 mag, in steps of 
0.1 mag. The number of synthetic reddened 
SEDs is around 2000. 

3. Definition of the merit functions: 

The best fits between observed and reddened 
synthetic SEDs in each aperture were found 
by using two merit functions. The first merit 
function compares the observed colors with 
the reddened SSP colors, and was defined as, 



x 



Wj (color° hh — color™~ dl ) 2 



(2) 



where the index i runs for 9 different col- 
ors, using as a reference the r band (FUV-r, 
NUV-r, u-r, g-r, r-i, r-z, r-J, r-H and r-K). 
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The variables color° hs and color™ dy are the 
observed and synthetic colors respectively, 
and <r° bs are the observational errors (instru- 
mental and external errors), the weight of 
each individual color is given by u>i. Some 
of the color indices are more sensitive to age 
and reddening than others, which could, in 
principle, be reflected by the corresponding 
weights. However, after several trials, we 
found that the best way to get a good fit 
(low~x 2 ) is to give the same weight for all the 
colors and introduce a separate merit func- 
tion to break the age-reddening degeneracy. 

The u-g color is chosen as a second merit 
function, given that the Balmer jump is 
the best discriminator between populations 
of different ages, for ages less than around 
500 Myr, and that it is least sensitive to ex- 
tinction. It is defined as: 



F=r [(u- 9 r s -(u-g) 



\mdl-\2 



(3) 



4. Selection of best-fit models: 

Due to the age-extinction degeneracy, dif- 
ferent combinations of age and color excess 
give us similar SEDs and therefore similar 
values of % 2 and p 2 . Instead of the common 
practice of retaining only the minimum x 2 
model, we retained 50 SSPs (out of the 2000 
reddened SSPs) with the lowest \ 2 ■ Of these 
50 SSPs, we then retained only those fits 
that are within 10% of the minimum value 
of (3 2 . Thus, for every aperture, we have a 
minimum of f and a maximum of 50 best-fit 
solutions. 

5. Monte-Carlo simulation of the observational 
errors: 

The error in each photometric band of an 
observed SED allows a set of synthetic SEDs 
to reproduce the observed SED with similar 
value of x 2 . In order to select all the syn- 
thetic SEDs that can fit the observed SED 
of an aperture, we generated 5000 new SEDs, 
where the flux in each photometric band was 
re-calculated by adding an error to the ob- 
served flux, where the error is randomly cho- 
sen such that their distribution is gaussian 
with a sigma equal to cr° bs for each band. We 
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Fig. 4. — The minimum \ 2 values obtained using 
Calzetti (solid circles) and Cardelli (open circles) 
extinction curves for all the selected apertures. 
Calzetti attenuation curve provides lower \ 2 va l~ 
ues, especially in the dusty central zone (apertures 
30-90). 



refer this set of new SEDs as Monte-Carlo 
SEDs or for brevity MC-SEDs. We repeat 
steps 3 and 4 with each of the 5000 MC- 
SEDs, and obtain the set of best-fit models 
(maximum of 50 models) for each one, re- 
sulting in a large set (several thousands) of 
best-fit combinations of ages and extinctions 
for each aperture. 

6. Determination of age and extinction: 

The age and extinction associated with each 
of the thousands of best-fit SSPs were ana- 
lyzed statistically using a gaussian distribu- 
tion function, to determine the most likely 
age and extinction of a region. Gaussian 
functions are fitted to the distributions of 
log(age) and E(B—V), with the center of the 
gaussians giving the most likely value, and 
sigma, an estimation of the error on the de- 
termined value. In a few cases where the dis- 
tribution has more than one peak, the gaus- 
sian function is fitted to the group having the 
maximum number of best-fits. The method 
is illustrated in Figure |3] for the brightest 
aperture in complex B. 

We found that the width of the distribution of 
ages and extinctions for a given region (defined by 
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Fig. 5. — Ages of the apertures obtained us- 
ing Calzetti (left) and Cardelli (right) attenuation 
curves. Both the curves give similar ages for disk, 
where as Cardelli curve gives systematically higher 
ages (~ 10 8 yr instead of ~ 10 6 yr ) for nuclear 
apertures. 



the <r of the best-fit gaussian) is within 15% of 
the derived value. The uncertainties were >30% 
if the observational errors in any of the bands is 
greater than 10%. As mentioned before, we dis- 
carded eight apertures based on that restriction. 
All these discarded apertures have FUV magni- 
tudes fainter than the 20a limiting magnitude. 

We adopted the Calzetti extinction curve in our 
analysis. Howev er, we have investiga ted the ef- 
fect of using the ICardelli et al.1 (| 19891 ) extinction 
curve instead. The x 2 values for the best-fit model 
for each aperture using the two different extinc- 
tion curves are compared in Figure |4j Both ex- 
tinction curves give similar values of % 2 for the 
disk apertures where the visual extinction is less 
than 2 mag. For the nuclear apertures Calzetti 
law gives systematically better fits. The ages de- 
rived for the disk apertures using the two curves 
do not differ much as can be seen in Figure [5] 
However we found that the ages derived using the 
Cardelli curve for majority of the nuclear regions 
are around 10 s years, which is clearly inconsistent 
with the estimation of ages using o ther methods 
(e.g. Forster Schreiber et al. I [20031) . Hence, we 
use the results obtained with the Calzetti curve 
throughout this paper. 



4. Analysis of Derived Ages and Extinc- 
tions in Named Complexes 

Table [T] lists the results obtained using the fit- 
ting method for apertures coincident with the 
high surface brightnes s com plexes defined by 



O'Conneh 1 k Manganol (|l978h . This table also 
shows spectroscopic ages and extinctions of star 
clusters located inside the corresponding aper- 
tures. Except for the complex B, the rest of the 
bright complexes have spatial extensions similar 
to or smaller than the aperture size, and hence the 
derived ages from aperture fluxes are the mean age 
of the complexes. The complex B encloses 8 aper- 
tures, with the aperture 27 being the brightest. 
Each of the complexes contain several compact 
star clusters. We used the high- resolution HST 
images to estimate the contribution of the clus- 
ters to the aperture fluxes, and found that the 
clusters contribute less than 1% to the disk aper- 
tures, whereas their contribution is less than 10% 
for the nuclear apertures (see §6). Hence, ages 
derived from the aperture fluxes do not represent 
those of the clusters, instead they are representa- 
tive of the field stellar populations of M82. The 
SEDs of the best-fit models for complexes A, B 
(aperture 27), C and D are shown in Figure [6j for 
selected apertures of complex B in Figure ; 7, and 
for complexes E, F, G and H in Figured] 

4.1. Nuclear complexes: A, C, D and E 

Complexes A, C, D, and E are located in the nu- 
clear star-forming region. This part of the galax y 



shows strong Ha emission ( Ohvama et al. I 2002) 



The presence of ionized gas poi nts to a population 
of stars of a few million years ( Rieke et al.lll98 
Satvapal et al.ll997tlForster Schreiber et al 11200 



More recently Konstantopoulos et al. ( 20091 ) found 
possible Wolf-Rayet features in the spectra of clus- 
ters located in these complexes, confirming that 
the clusters of the nuclear complexes are young. 

Complex A is the core of the nuclear star-burst 
and is the brightest complex. This is the most 
studied complex of M82. Surrounding the complex 
A are complexes C, D and E. Each one of these 
complexes are resol ved into several tens of CSCs 
in the HST images (|0'Connell et al. Ill995l ). Our 



aperture 57 encloses the bright core of the com- 
plex A, for which we obtained an age of ~ 2 Myr 
(log age — 6.31 ± 0.05) and A v = 5.7 mag. 
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Fig. 6. — The observed SEDs for complexes A, B, C and D (solid circles) are compared with the best-fit 
model SED (open circles) and the corresponding full resolution model spectrum (solid line). The derived 
age and extinction are given in each plot . 



For complexes C and E, we derived ages of 
around 5 and 7 Myr (log age = 6.66 ± 0.03 
and log age — 6.85 ± 0.02, respectively) with 
A v values of ^5 mag. These derived param- 
eters are consistent with the previous results 



(|F5rster Schreiber et al. 1 120031 ) . For complex D, 



which is directly to the north of complex A, we de- 
rived an age of ~95 Myr (log age — 7.97 ± 0.07). 
This age is clearly larger than the age of the nu- 
clear starburst and probably this belongs to the 
disk and it is seen projected on top of the nucleus 
of the galaxy. 

4.2. Complex B 

Complex B is the most studied complex out- 
side the nuclear zone. It lies around a kilopar- 
sec to the north-east of the nucleus. In the blue 
band images this complex is even brighter than 
the nucleus. Complex B spans over a large area of 
the disk of M82, and is the largest of the named 



complexes. Sometimes, it is subdivided into two 
parts, the Bl containing the bright zones to the 
far east and B2 containing relative ly fainter zones 
to the west ( de Griis et al.l 200lh. U sing long- 
slit spectroscopy, Mavva et al. ( 2006f ) obtained 
an age of 500 Myr for this complex. There ar e 



40 CSCs in this complex (jMavva et al. 1120081 ) 



Recent spectroscopic studies have obtained ages 



for some of these clusters. ISmith et al.l (|2007h 
found that stellar clusters located in this area 
cover an ag e interval of 80-270 M y r, pea king at 



150 Myr. iKonstantopoulos et al. (|2009h deter- 



mined that star clusters located in this complex 
are 80-200 Myr old. 

In our analysis, complex B is covered by 8 aper- 
tures. We have listed these apertures in Table[T]as 
function of the distance to the center of the galaxy. 
Figure [7] displays the best-fit SEDs for selected 
apertures located inside this complex. Apertures 
27 and 29, which lie to the extreme east contain 
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Fig. 7. — Similar to Fig. [6j but for some apertures located inside complex B. 



the brightest part of this complex, for which we 
derived ages of -400 Myr (log age = 8.60 ± 0.06, 
and log age — 8.59 ± 0.06, respectively). The 
age derived for the aperture 50 - which lies on 
the extreme west of this complex - is —180 Myr 
(log age — 8.25 ± 0.07). We note that the ages of 
the western apertures are systematically younger 
reaching values of around 200 Myr at the extreme 
west (see Tabled]). 

4.3. Complexes F, G and H 

Complexes F and H lie in the transition zone be- 
tween the nucleus and the disk. Each of these com- 
plexes is associated with a very luminous and mas- 
sive stellar cluster. Complex F is located on the 
south-west of the nucleus (aperture 72), whereas 
complex H is at the north-east of the nucleus 
(aperture 52). Complex G (aperture 80) is located 
in the disk on the west side. This latter complex 
was not included in previous spectroscopic studies 
of M82 probably because it is very diffuse. 



Smith fc Gallagherl (|200ll ) performed a spectro- 
scopic study of the stellar cluster of complex F. 
They found that this cluster is around 50 Myr old. 



Bastian et al.l (|2007l ) found substructures within 



the cluster, and in the surrounding region, which 
was interpreted as due to differential extinction 
across the face of this complex, with some re- 
gions having A, > 2.5 mag. They estimate d an 
age of 60 Myr. Konstantopoulos et al.l ( 20091 ) an- 
alyzed two clusters in this area (slits 43.2 and 
51(F)). They concluded that the age of cluster F 
is 40-80 Myr (in accor d ance w ith r esult obtained 



by | Smith fc Gallagherl (|200ll ) and iBastian et al 
j2007j)), whereas cluster 43.2 is 100-320 Myr old. 



Complex H is dominated by one of the most 
massive star cluster of M82 (~ 2.6 x 10 5 M , 
Mavva et al. ( 20081) 1 which is also known as clus- 



ter B2 -1. Using ffffTimages. lKonstantopoulos et al 



(|20081 ) found a peculiar extinction pattern asso- 
ciated to this cluster, complicati ng its age determi- 
nation . In a more recent analvsis lKonstantopoulos et al 
(2009) determined that the cluster of complex H 
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Fig. 8. — Similar to Fig. |6l but for complexes E, F, G and H. 



is 180-200 Myr old. 

Figure [8] displays the best-fit models for these 
complexes. For complex F we found a stellar pop- 
ulation of around 180 Myr (log age = 8.25±0.07). 
For the population of complex G we determined an 
age of ~ 350 Myr. For the complex H, we obtained 
an age of ~ 75 Myr (log age = 7.87 ±0.05), which 
is younger than the star cluster B2-1. 

5. Spatially Resolved Star Formation His- 
tory and Extinction 

In this section, we use the age and extinction 
of all the apertures to infer the spatially resolved 
star formation history in the disk of M82 in scales 
of - 200 pc. 

Figure [5] shows the distribution of age as a func- 
tion of the distance to the center of the galaxy. 
We distinguish those apertures belonging to com- 
plexes A to H by denoting them with filled cir- 
cles. Two distinct epochs of star formation can be 
easily inferred from this plot: one younger than 



10 Myr, exclusively located in the central 0.5 kpc 
radius, and the other older than 100 Myr, dis- 
tributed along the whole disk of the galaxy. In 
Figure [TU1 we display the age of the apertures as 
a function of the distance using linear scale. A 
trend can be seen in this plot in such a way that 
the ages decrease from around 450 Myr at ^3 kpc 
to 300 Myr at 0.5 kpc on the eastern side (positive 
distances) , and from 400 Myr at 3 kpc to 300 Myr 
at 0.5 kpc on the western side. In the transition 
zones between the disk and the nucleus, the ages 
decrease rapidly from ~ 300 Myr to < 100 Myr. 
The observed trend may correspond to a smooth 
age gradient in the disk. However, the uncertain- 
ties in the derived ages in this work, do not permit 
a clear answer on this. 

Another important thing to note in Figure [TOl is 
the similarity in ages of the high surface brightness 
named complexes and the surrounding disk. This 
indicates that the stellar populations in the field, 
as well as in the bright complexes, were formed in 
the same event. Implications of these results will 
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Age of the Stellar population and extinction found in apertures related with complexes 
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References. — (1) Smith et al. 2006; (2) Bastian et al. 2007; (3) Smith et al 2007. 
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Fig. 9. — Derived ages of the dominant stellar pop- 
ulation within apertures of 5" radius as a func- 
tion of the distance to the galactic center (posi- 
tive distances correspond to the north-east direc- 
tion). Circles mark the results for the selected 
apertures, with th e filled circles marking t h e com - 
plexes defined by lO'Connell fc Manganol ll 19781) 



Ages of the clusters from iKonstantopoulos et al 



( 2009ft are denoted by triangles 



be discussed in §7. 

Konstantopoulos et al. ( 2009ft used multi-object 



spectroscopy to obtain ages and extinctions for 49 
relatively bright clusters, covering up to a radial 
distance of 3 kpc. Their ages are displayed in Fig- 
ure IH1 using filled triangles. They found a mean age 
of 150 Myr for disk clusters, with no noticeable 
age gradient. Hence, these clusters seem to be 
systematically younger than the field populations 
as traced by our apertures. This difference makes 
the clusters slightly bluer than the surro unding 



disk as was noted by iMavva et al~l (|2008ft . The 



implications of these differences are discussed in 
§7- 

Several studies have found an extinction gra- 
dient in the central kiloparsec of M82 with 
the extinction inc reasing towards the nucleus 
( Waller et al.l Il992i ). Most of these studies are 



based on nebular line ratios and hence sample 
only the central part of M82. Our method al- 
lows mapping of the extinction in the entire disk. 
In Figure II H we show the radial distribution 
of extinction. A gradient of extinction can be 
seen in the disk on both sides of the nucleus, 
with the A v increasing smoothly from 1.5 mag 




2000 1000 -1000 -2000 
Distance (pc) 

Fig. 10. — As in previous Figure, but with the 
vertical axis on a linear scale to enable better vi- 
sualization of a slightly outside-in age gradient in 
the disk of M82. 



at 3 kpc to 2.8 mag at 0.5 kpc. For the nuclear 
apertures, we found extinction ranging between 
3-6 mag, which agree s well with the values ob- 
tained bv lWaller et "aH (Il992l ) using nebular lines. 
On the other hand, the star clusters studied by 
Konstantopoulos et al. I (|2009ft present extinction 
values that are systematically ~ 1 mag lower than 
our values for selected apertures (see Tabled]). 



5.1. The underlying old disk 

It is important to note that we could reproduce 
the SEDs of all the disk apertures using a single 
stellar population of age ~ 400 Myr, without the 
need to add an SSP of several gigayears of an un- 
derlying old disk. These results are consistent with 
the disk-wide star-fo rmation episode proposed by 
Mavva et al~l (|2006ft . in which more than 90% of 
stellar mass in the disk was formed. The derived 
ap erture age s are co nsistent also with that inferred 
by iDavidge" ( 2008ft . a nd are marginally yo unger 
than that proposed by IMavva et al. ] (12006ft . 

In this section, we check whether our multi- 
band SEDs are consistent with less than 10% of 
mass in old underlying disk stars. In order to 
do this, we performed new fits between the ob- 
served SEDs and new models. The new models 
were calculated by adding an SSP of 4 Gyr old to 
the best-fit SSP discussed in previous sections. We 
assumed that the old component shares the same 
reddening as the best-fit SSP. We carried out two 
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Fig. 11. — The visual extinction within apertures 
of 5" radius as a function of the distance to the 
galactic center (positive distances correspond to 
the north-east direction). The central part of 
the galaxy is more extincted, Ay > 3.5 mag. 
Circles describe our results, filled circles are re- 
lated to named complexes, and triangles represent 
Ay in clusters obtained bv lKonstantopoulos et al 
(|2009h . 



tests: one using solar metallicity for the old com- 
ponent and a second test using a metal poor model 
(Z= 0.008) for the old component. The mass of 
the old population (old-pop) model is varied from 
a minimum of 0% in mass of the best-fit SSP in 
steps of 1% until the flux of the new model ex- 
ceeds the observed flux plus the estimated error in 
the K band. We have chosen the K-band, because 
the old population is expected to have a maximum 
contribution in this band. The analysis performed 
for an area of complex B (aperture 27) using solar 
metallicity is illustrated in Figure Q21 The upper 
plot displays the observed SED (dots), the SED 
of the best-fit model (solid line fitting the dots) 
and the SED of the old population (solid line) 
when it has a 9% in mass of the best-fit model. 
The lower plot shows the comparison between the 
flux of the old-pop model (with different percent- 
age of mass) and the flux of the best-fit model 
(i?atio = F old _ Pop /F Best _ Fit ). It can be inferred 
that the old population with solar metallicity at 
the most contains 9% of the stellar mass of the 
best-fit population inside aperture 27. In the case 
of a metal-poor model, the upper limit in mass is 
11%. More massive old stellar population would 



have increased the K-band flux to values above 
that observed. The upper mass limit of the old 
population obtained for aperture 27 is typical for 
the rest of the disk apertures, as is illustrated in 
the histogram of Figure (T3J The mean upper mass 
limit is 9% and 11% depending on whether the 
old population is assumed to have solar or sub- 
solar metallicity, respectively. Thus, the observed 
multi-band SEDs of individual aper tures are con- 



sistent with the result obtained by iMavva et al. 



( 2006( ) where the underlying old disk, if present, is 
very insignificant in M82. 
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Fig. 12. — The observed SED (circles) of an area 
of complex B (aperture 27) is plotted in the up- 
per panel, as well as the SED of the best-fit SSP 
(solid line fitting the dots) and the SED of the 
old population model (solid line) with 9% in mass 
of the best-fit model. The lower panel shows the 
ratio {F id- Pop / F Bes t- Fit) for different fractional 
masses (indicated as percentage for each curve) of 
the old population. For the aperture plotted, a 
mass>9% would have been noticed in the K-band 
photometry. 



6. Cluster Formation Efficiency 

From the analysis of the previous sections, it is 
clear that the observed aperture fluxes of both the 
disk and nuclear regions of M82 are dominated by 
a single generation of stars modeled as a simple 
stellar population. M82 contains a rich popula- 
tion of CSCs, with each CSC containin g a single 

l2008h . So, it is 



generation of stars ([Mavva et al. 



interesting to investigate whether the CSCs and 
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Fig. 13. — The histogram shows the limiting upper 
masses of old population that can be added to the 
best-fit SSP of apertures located in the disk of 
M82. The mean value is 9% for solar metallicity 
models and 11% for metal-poor models. 

the field stellar populations formed from the same 
episode. 

In Figure G3 we have compared the ages derived 
using our aperture fluxes to the ages of stellar clus- 
ters located within the apertures. The disk clus- 
ters have a mean age of 150 Myr, whereas most 
of the field populations surrounding these clusters 
have ages between 300-400 Myr. Thus, clusters 
seem to be systematically younger than the sur- 
rounding field populations by around 200 Myr. 
In the disk-wide st arburst model proposed by 
Mav va et al. (200G), the star formation episode 
in the disk lasted for around 300 Myr. Thus, in 
spite of the systematic differences in mean ages, it 
is still possible that both the clusters and the field 
stars formed from the same episode. On the other 
hand, in the nuclear region, both clusters and stars 
in their vicinity have similar ages, both popula- 
tions being part of the on-going nuclear starburst. 

Another question of interest is to find out 
whether clusters dominate the observed fluxes 
at spatial scales of 180 pc. In order to answer 
this question, we compared the total fluxes of all 
clusters located physically within an aperture of 
5" radius with the integrated fluxes within that 
aperture. Typically, nuclear apertures contain 
more than 10 compact clusters with the aper- 
ture 57 (complex A) containing as much as 76, 
whereas disk apertures contain less than 10 clus- 
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Fig. 14. — The contribution to the total flux by 
the stellar clusters as a function of the distance to 
the galactic center (each point corresponds to an 
aperture ). The upper plot shows that the mean 
contribution in the nucleus is ~ 7%, whereas in the 
disk it is around 0.7% (dashed lines). The lower 
plot shows the number of star clusters inside each 
aperture (as in previous figures, positive distances 
correspond with the NE). 



ters and several apertures do not have any cata- 
logued cluster inside. We calculated the fluxes of 
the field population using the apertures defined 
in Figure 1 on the F555 W image obtained from 



EST /ACS observations (|Mutchler et al. I [2007) 



The F555W band c luster fluxes were taken from 
iMavva et al~l (|2008l) . 

Figure [14] displays in the upper plot the com- 
parison between the observed flux and the fluxes 
emitted by stellar clusters in the selected aper- 
tures. The lower plot shows the number of star 
clusters inside each aperture. Clusters contribute 
clearly more to the aperture fluxes in the nucleus 
as compared to the disk, 7% vs. 0.7%. These re- 
sults show that field stars, not the clusters, are 
emitting most of the observed flux at scales of 
around 180 pc, even in areas where there is strong 
star formation. Similar results are obtained with 
the F435W and F814W images from EST /ACS. 
There are apertures located inside the 500 pc ra- 
dius with ratios below ~2%. These apertures seem 
to belong to the disk far away from the center of 
the galaxy but due to the high inclination of the 
galaxy are seen projected on to the nuclear region. 
The derived ages and extinctions of these regions 
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confirm that interpretation. 

The flux fraction obtained above is directly re- 
lated to the fraction of stellar mass in compact 
clusters as compared to the mass in all field stars, 
provided the clusters and field stars are of the 
same age. For the nuclear region this assump- 
tion holds, and hence on an average, 7% of the 
stars in the nuclear region of M82 form in clus- 
ters. As these nuclear populations are younger 
than 10 Myr, the observed fraction corresponds 
to the cluster formation efficiency. Curiously, this 
estimate compares very well with the clu ster for- 
matio n efficiency of 8 ± 3% determined by lBastianl 
(2008) in star-forming galaxies. 

On the other hand, the disk clusters are younger 
than the surrounding field stars. Given that the 
fluxes of stellar population decrease with time due 
to passive evolution, the observed flux fraction of 
0.7% is an upper limit to the fraction of mass in 
cluster stars. 



Meurer et al. 



(|1995[ ). working with the HST 
UV images of nine starburst galaxies with several 
bright star clusters, found that 20% of the total 
UV luminosity is emitted by the clusters. It is 
most likely that the typical age of star formation 
in his sample galaxies is only around 10 Myr. Our 
values for the nucleus of M82 reaches 20% in one 
case with the rest havin g valu es below those ob- 



tained bv I Meurer et al. I (J1995). So it seems that 



CSCs contribute < 20% to the flux in young star 
forming regions decreasing to less than one percent 
when they are several hundreds of million years 
old. 

7. Discussion and Implications 

We start our discussion addressing the na- 
ture of complex B, which has drawn a lot of 
attention in recent years (|de Griis et i 



Smith et al.l 120071 : iKonstantopoulos et al 



2001 



2009b|) 



Often this complex is refer red to as fossil star- 



burst (jde Griis et al.l 120011) . There was also 



the suggestion that the brightness of this com- 
plex is due to windows in th e dust distribution 
I Konstantopoulos et al. 20081 ). Results obtained 



with our fitting method show that ages and ex- 
tinctions of the bright named complexes are simi- 
lar to the surrounding relatively fainter areas (see 
Figures [9land [TTj) . suggesting that the high bright- 
ness of the named complexes is due to their higher 



star formation rate with respect to the surround- 
ing disk, and not because of lower extinction or 
younger age. Particularly, apertures belonging to 
complex B do not show significant lower extinct ion 



as often claimed ([Konstantopoulos et al. 2008) 



Ages of CSCs obtained us ing spectroscopy by 
Konstantopoulos et al. ( 20091 ) are in the range of 



100-300 Myr with the mean age of 150 Myr. Thus 
clusters seem to be systematically younger than 
the surrounding disk stars. However, it is impor- 
tant to point o ut that only around 50 of the 653 
known clusters (|Mavva et al. 1120081 ) were selected 
for spectroscopic observations and hence the de- 
rived ages may not be representative of the entire 
cluster population. Given that the cluster lumi- 
nosity decreases with age, the brightest clusters 
are expected to be the youngest. Hence it is not 
strange that the derived ages of the clusters are 
systematically younger. Complex H, which lies in 
the transition zone between the nucleus and disk, 
offers an interesting case, where the compact clus- 
ter is older than the field population surrounding 

it. 

Mavva et al. ( 20061 ) found that the entire disk 



participated in a star-forming episode that started 
following the interaction with M81 around 500 
Myr ago. Ages of the dominant population found 
in our study are consistent with that idea. How- 
ever, we find in the disk of M82 a systematic trend 
for the external zones to be marginally older as 
compared to the inner zones. Considering that our 
ages are luminosity weighted, this trend suggests 
either that the star-forming episode was not simul- 
taneous in the entire disk but it is systematically 
younger in the inner regions, or that the entire 
disk started forming stars at the same time, but 
the star-forming episode stopped systematically at 
earlier epochs in the outer regions. Study of re- 
solved stellar populations bv lDavidge" (2008) sup- 
ports the latter idea over a much extended disk. 

We commented that the mean SEDs of central, 
north-east and south-west parts of the disk dis- 
play differences (see Fig. [2]). The age and extinc- 
tion distributions along the disk of the galaxy can 
explain these differences. The visual absorption 
found on each side of the galaxy is similar, the 
mean Ay of the NE is 1.6 mag and the mean Ay 
of the SW is 1.8 mag. However, the population in 
the south-west is slightly younger: the oldest pop- 
ulation of the SW is around 400 Myr old whereas 
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in the NE side the oldest population is ^450 Myr 
old (see Figs. [TUl and ITTj) . This behavior explains 
the bluer appearance of mean SED of the south- 
west part. On the other hand, we have pointed out 
that the mean SED of the central part is redder 
than the mean SEDs of the NE and SW regions. 
In spite of the very young ages found in the cen- 
tral parts of the galaxy, the high optical extinction 
(Ay around 5 mag, see Fig. found in this part 
produces a mean SED which looks very reddened. 

There is very little evidence for the existence of 
stars formed previous to the known star-forming 
episode. In the star-formation mod el proposed to 



explai n the properties of the disk, iMayya et al. 



( 20061 ) found a maximum contribution of around 
10% of the total mass in stars older than 1 Gyr. 
In our study, we have analyzed this fraction in 
every aperture. We found that old stars do not 
contribute more than 11% in majority of the disk 
apertures. Thus, we reiterate that the disk of M82 
formed almost all of its stars after the interaction 
with M81. 

Young compact clusters contribute ~7% to the 
nuclear aperture fluxes, with the corresponding 
contribution being ^0.7% for the relatively older 
disk apertures. Given that the clusters and the 
surrounding regions share similar age and extinc- 
tion, the above fraction is an indication of the ef- 
ficiency of cluster formation and its dissolution. 
Hence for young nuclear clusters, the efficiency on 
an average is around 7%. The observed relatively 
lower fractions for the disk regions may indicate 
that on the time scale of a few hundred million 
years many clusters are already dissolved. 

8. Summary 

In this work, we analyzed the star formation 
history of M82 over spatial scales of ~ 180 pc 
using the broad-band SEDs covering the FUV to 
NIR bands. We constructed 117 SEDs, where each 
point of the SED is obtained by carrying out pho- 
tometry in apertures of 5" radius over the entire 
face of M82 within 3 kpc radius. These SEDs were 
fitted with synthetic SSPs in wide range of red- 
dening and ages. With this technique, we con- 
firmed the following two facts that were already 
known: (1) the stellar populations in the nuclear 
zone are less than 10 Myr old and are heavily red- 
dened (A v =2-6 mag), (2) the stellar populations 



in the disk have ages between 100-450 Myr. 

The use of spatially resolved SEDs has al- 
lowed us to establish, for the first time, that the 
interaction-driven disk star formation was not re- 
stricted to a few bright zones, but occurred every- 
where in the disk over spatial scales as small as 
180 pc. Each of these SEDs is consistent with the 
idea that more than 90% of the stellar mass of the 
disk is formed as part of the disk-wide star forma- 
tion episode, and containing only a small fraction 
of stars older than 1 Gyr. 

We found that the bright named complexes 
have similar ages and reddening as that of the 
surrounding low-surface brightness parts of the 
galaxy: complexes within the starburst region 
(complexes A, C, E) have ages less than 10 Myr 
and the complex B in the disk has age 200- 
400 Myr, and the complexes in the intermedi- 
ate zone (F and H) having ages between 50- 
200 Myr. This implies that the named complexes 
are brighter just because of higher amount of star 
formation, not because they are younger or have 
lesser reddening. Ages of compact stellar clus- 
ters in the disk derived using spectroscopy (~ 
150 Myr) are systematically smaller than the mean 
age of the field stellar populations (~ 350 Myr). 
However, the difference in the age is of the same 
order as the duration of the disk-wide star forma- 
tion episode, and hence both populations could 
belong to the same episode. 

In the nuclear regions, around 7% of all the stel- 
lar mass resides in the cluster stars, whereas that 
fraction is less than 1% in the disk. Given that 
the disk populations are systematically older than 
the nuclear populations, the difference in the frac- 
tional mass in cluster stars may be due to disrup- 
tion of clusters as they evolve. 

The results obtained in this study are encour- 
aging from the point of view of extending the tech- 
nique adopted in this work to other nearby galax- 
ies. Thus star formation history over a few hun- 
dreds of parsec scale can be derived by making use 
of the already available multi-band archival images 
such as GALEX, SDSS, and 2MASS. 
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